A Majorana smoking gun for the superconductor-semiconductor hybrid topological 

system 



S. Das Sarma 1 , Jay D. Sau 2 , and Tudor D. Stanescu 3 
1 Condensed Matter Theory Center, Department of Physics, 
University of Maryland, College Park, Maryland 20742-4111, USA 
2 Department of Physics, Harvard University, Cambridge, Massachusetts 02138, USA 
3 Department of Physics, West Virginia University, Morgantown, WV 26506, USA 

Recent observations of a zero bias conductance peak in tunneling transport measurements in 
superconductor-semiconductor nanowire devices provide evidence for the predicted zero-energy 
Majorana modes, but not the conclusive proof for their existence. We establish that direct ob- 
servation of a splitting of the zero bias conductance peak can serve as the smoking gun evidence 
for the existence of the Majorana mode. We show that the splitting has an oscillatory dependence 
on the Zeeman field (chemical potential) at fixed chemical potential (Zeeman field). By contrast, 
when the density is constant rather than the chemical potential - the likely situation in the current 
experimental set-ups - the splitting oscillations are generically suppressed. Our theory predicts the 
conditions under which the splitting oscillations can serve as the smoking gun for the experimental 
confirmation of the elusive Majorana mode. 

PACS numbers: 



The recent experimental report pQ providing direct 
observational evidence for the possible existence of the 
predicted [2}{5] zero-energy Majorana quasiparticle in 
superconductor-semiconductor nanowire hybrid struc- 
tures in the presence of spin-orbit coupling and Zee- 
man splitting has tremendously excited the whole physics 
community [6}{9]. Yet, in spite of several subsequent re- 
ports [TQHT2] having validated the original data of Ref. 
[TJ this experiment has also raised many questions. Most 
of these questions arise from a critical comparison be- 
tween the experimental data pQ and the original theoret- 
ical predictions [2HS] , leading to the inevitable conclusion 
that there are some significant discrepancies between ex- 
periment and theory. For example, the key experimental 
observation is the development of a robust subgap zero 
bias conductance peak (ZBCP) in the tunneling differen- 
tial conductance of the nanowire in the presence of an ex- 
ternal magnetic field applied along the wire, as predicted 
theoretically [3] and as expected for a Majorana zero en- 
ergy mode in a topological superconductor [l3j[T4]. How- 
ever, the actual magnitude of the ZBCP (~ 0.1e 2 /h) is 
more than an order of magnitude smaller than the pre- 
dicted ideal quantized value (2e 2 /h). In addition, the 
Major ana-induced ZBCP should only appear beyond a 
magnetic field-driven topological quantum phase transi- 
tion (TQPT) characterized by the closing of the super- 
conducting (SC) gap, yet there is no apparent signature 
of gap closing in the measured tunneling current. Al- 
though recent theoretical works [T5HT7] provide reason- 
able explanations for some of these discrepancies, other 
recent papers emphasize that a ZBCP could arise in the 
system in the absence of Majorana bound states, due 
to more mundane mechanisms involving strong disorder 
[T8H2Q] , smooth end confinement [21], or Kondo physics 
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Given this fluid and confusing nature of the subject 
matter, with publications arguing in favor of or against 
the Majorana interpretation of the experimental observa- 
tion in Ref. 1 appearing almost weekly, it is of paramount 
importance to conceive of hallmark experimental signa- 
tures for the Majorana quasiparticle. It was already em- 
phasized in the original theoretical predictions [3 , 4 that 
the observation of a ZBCP at finite magnetic field is only 
a necessary condition for the existence of Majorana quasi- 
particles. The sufficient condition to validate their exis- 
tence must be some type of interference measurement, 
such as the fractional Josephson effect [23] [24] manifest- 
ing a 47r periodicity in an ac Josephson measurement. 
While such a measurement, or the direct observation 
of non-Abelian Majorana interference [25H27] . will cer- 
tainly be necessary down the line to absolutely validate 
the existence of localized non-Abelian Majorana modes, 
the high level of complexity and difficulty of this type 
of measurements make it unlikely that they will be suc- 
cessful in the near future. Therefore, it is desirable that 
a simpler experimental smoking gun for the Majorana, 
something with a difficulty level comparable with to the 
existing ZBCP experiments, be proposed and carried out 
long before the rather challenging fractional Josephson 
measurement and the interference experiments could be 
performed. In the current work, we propose a smoking 
gun Majorana measurement that could not only be car- 
ried out right now, but in fact it is conceivable that the 
necessary experimental data could already be hidden in 
the reported ZBCP measurements. 

Our key observations in this context are: i) nanowire 
Majorana modes always come in pairs [23] localized at 
the two ends of the wire, and ii) the Majorana mode is 
a pure zero-energy mode only when the wire is infinitely 
long. For any realistic finite-length wire, the two end 
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Majorana wavefunctions overlap and the resulting hy- 
bridization leads to a splitting of the zero- mode [28 . This 
hybridization-induced energy-splitting increases with ei- 
ther decreasing wire length or increasing coherence length 
and is characterized by an oscillatory behavior deter- 
mined by the Fermi wavevector of the top occupied band 
[28] . This is in contrast to other sources of zero-bias 
conductance peaks such as Kondo resonances [22] where 
the Zeeman splitting is expected to produce a monoton- 
ically increasing splitting. Since both the effective SC 
coherence length and the effective Fermi wavevector can 
be tuned by varying the external magnetic field or the 
chemical potential of the system, a direct observation of 
the oscillatory energy splitting with increasing magnetic 
field (which increases the coherence length by increasing 
the effective Fermi-velocity and suppressing the SC gap) 
should be a definitive smoking gun for the Majorana ex- 
istence. 

In real systems, there are two main challenges that 
have to be addressed: the finite energy resolution and 
the limited experimental ability to control the chemical 
potential. Various broadening mechanisms (e.g., temper- 
ature, inelastic scattering, quasiparticle poisoning, disor- 
der, finite tunnel barrier) may mask the underlying Ma- 
jorana splitting oscillations as two closely-spaced split 
peaks near zero bias may merge into a single broad zero- 
bias peak. Nonetheless, the width of the Majorana peak 
should show a modulation with varying Zeeman field or 
chemical potential, even if the finite resolution masks the 
splitting itself. In the current work we make detailed the- 
oretical predictions about how the Majorana ZBCP split- 
ting should depend on the relevant parameters. These 
predictions are experimentally verifiable and can serve as 
a clear smoking gun for the validation (or invalidation) of 
the observation of Majorana modes in hybrid nanostruc- 
tures. Since four distinct experimental groups [TJ HQUT2] 
have already reported detailed results for the ZBCP de- 
pendence on the applied Zeeman field, it appears that our 
proposed Majorana smoking gun measurement should be 
feasible in the very near future. 

The currently existing experimental results do not 
manifest any clear signature of Majorana oscillatory 
splitting, although the issue is by no means definitively 
sorted out since the higher field data are rather sparse 
and have not yet been carefully analyzed. However, 
more importantly, our work shows that one must con- 
ceptually distinguish between the field-dependent Majo- 
rana splitting at constant chemical potential versus con- 
stant carrier density. All discussions of Majorana split- 
ting in the literature have so far considered only the 
constant chemical potential situation, where both the 
average splitting and the oscillation amplitude increase 
monotonically with increasing Zeeman splitting, yet ac- 
tual experiments, because of repulsive interactions are 
mostly carried out at constant density. We show that 
the constant density condition is qualitatively different 
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FIG. 1: (Color online) Oscillatory splitting of the Majo- 
rana mode as a function of Vz (for fi = /is, top panel) 
and Afi = \i — \±z (for Vz = 0.75 meV, bottom). The 
parameters used in the calculation are: A = 0.25 meV, 
a 2 ra*/2/i 2 = 50/iey, and m e ff = 0.015ra e . This behavior 
is independent on the number of occupied bands. Finite en- 
ergy resolution (SE = 10/ieV, yellow band) may prevent the 
observation of oscillatory slitting in long wires. 

and may lead to a strong suppression of the oscillatory 
splitting, although generically some splitting is always 
present. Moreover, the oscillatory splitting is always 
present at constant Zeeman field as a function of chem- 
ical potential (or particle density). Thus, experiments 
searching for the oscillatory splitting (or the associated 
modulations of the ZBCP width) in gated samples should 
focus on the ability to tune the density and/or chem- 
ical potential. We know of no other proposed mecha- 
nism in superconductor-semiconductor hybrid structures 
that could lead to a ZBCP manifesting oscillatory split- 
ting, except for the Major ana-induced zero bias peak 
predicted to occur in the topological SC phase. 

Finite SC nanowires in the topological phase have been 
proposed to support a pair of near zero-energy Majorana 
fermions (MFs), one at each end of the wire. The energy 
degeneracy of these zero modes is removed [28] by an 
energy splitting 



AE ~ Re[V](L/2) Jtf r (Z,/2)], 



(1) 



where ^ j7 .(L/2) are the wave- functions of the MF bound- 
states near the left and right ends respectively and J is 
the splitting operator whose explicit form is given in the 
Supplementary Material. The spatial dependence of the 
MF wave-functions, which can be approximated as (see 
Supplementary Material for details) 

9i(x) oce- x/ Ze ±ikF ^f x , (2) 

for x ^> £ where £ is the effective coherence length and 
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kp, e ff is the effective fermi- wave- vector associated with 
the zero- mode solution. The parameter kp, e ff an d £ de- 
pend on the microscopic parameters such as the Zeeman 
splitting Vz and the chemical /i, which can be tuned ex- 
ternally by changing the applied in-plane magnetic field 
or gate voltage. The wave-function & r (x) of the MF at 
the right end behaves in a qualitatively similar way, so 
that for L ^> £ the MF splitting wave-function has the 
approximate form 

e~ 2L/ t 

AEzih kF,eff m ^ cos (k Fie f fL) , (3) 

where m is the effective electron mass in the nanowire. 
The cos (kp^ffL) factor should lead to an oscillation of 
the energy splitting AE as a function of the separation L 
between the MF. In the nanowire case, it is more natural 
to tune the Zeeman potential Vz or the gate voltage /i. 
Since kp, e ff depends on Vz and /i, for sufficiently long 
L, one can expect the cos {fcp je //(Vz, aO^} factor to lead 
to oscillations in AE as a function of both Vz and /i. 
We note that the e~ 2L ^ factor leads to the well-known 
exponential topological protection of the zero-energy MF, 
which, however, only applies in the very long- wire limit 
of L ^> £ where Eq. (3) is valid. 

To confirm these expectations, we study the spectrum 
of a spin-orbit coupled semiconductor nanowire in prox- 
imity to a superconductor and placed in a magnetic 
field parallel to the wire using the numerical diagonal- 
ization of the Bogoliubov-de Gennes (BdG) Hamiltonian 
[29] . Since we are considering an experimentally relevant 
quasi-lD nanowire (with diameter ~ 80nm), our system 
can support several sub-bands with dispersion minima 
near the spin-degenerate energies E n = E n (k = 0,V Z = 
0). As the energy separation between the sub-bands is 
larger than the SC gap, the top-most occupied band n top , 
with effective chemical potential A/i = /i — /i ntop , where 
/i nt0 p = ^n t0 p > dominates the low energy part of the spec- 
trum, including the properties of the MF bound states. 
The splitting of the MFs for a nanowire with n top = 3 as 
function of Vz and A/i = /i — /i3 is shown in Fig.[l]for two 
different realistic wire lengths. As expected from Eq. |3j 
the energy splitting AE is oscillatory and the oscillation 
period becomes shorter in longer nanowires, while the 
amplitude of the oscillations decreases. The exponential 
decrease of the oscillation amplitude is controlled by the 
factor e~ L l^ in Eq. [3J which reflects the exponential de- 
cay of the MF wave-function away from the wire-ends 
as in Eq. [2] This suppression of the energy splitting in 
long wires depends critically on the actual magnitude of 
the coherence length £. As seen in Fig. gb), the calcu- 
lated coherence length £ reaches a minimum for a value 
of Vz near the critical Zeeman field Vz ^ Vz, c , then in- 
creases monotonically as a function of Zeeman field. The 
increase in the coherence length £ can be attributed to a 
general decrease in the gap (shown in Fig. |2^a)) and a si- 
multaneous increase of the fermi velocity with increasing 
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FIG. 2: (Color online) Top: Dependence of the SC gap on 
the Zeeman field in the topological SC phase (yellow) for two 
different values of the chemical potential. Bottom: Decay 
length of the MFs as a function of Zeeman splitting. Note that 
the minimum of the coherence length occurs at Vz > Vz,c, 
right above the critical field, and does not coincide with the 
maximum of the energy gap. 

field. In addition, we find a weak increase of the coher- 
ence length with A/i. This dependence of the coherence 
length on Vz and A/i is ultimately responsible for the 
increase of the oscillation amplitude shown in Fig. [I] 

What are the main challenges facing the experimental 
detection of this unique signature of the Majorana mode, 
the oscillatory splitting of the ZBCP? First, the broaden- 
ing of the ZBCP implies finite energy resolution, hence in 
long enough wires, where the splitting is small, the oscil- 
latory splitting can be resolved only above a certain min- 
imum value of the Zeeman field. For example, in the case 
of a ZBCP with a broadening of 10 /ieV (shown as the 
yellow band in Fig. [TJ, the oscillations of the Majorana 
mode become visible in a ~ 1 /im wire for Vz > 0.5 meV. 
The required Zeeman field for a longer wire (L ~ 2 jam) 
is significantly larger (see Fig.[T]). However, reducing too 
much the length of the wire (to enhance the amplitude 
of the splitting) results in oscillations with very long pe- 
riod, which might again prevent the observation of the 
oscillatory behavior. Second, the oscillatory splitting of 
the Majorana mode as function of V z is a necessary fea- 
ture in a system with fixed chemical potential. However, 
this condition may not be satisfied in the current exper- 
iments. We emphasize that, while there appears to be 
a splitting of the ZBCP in the recent measurements pQ, 
the are no obvious oscillations of this splitting. Here we 
provide a mechanism which may lead to a considerable 
suppression of the splitting oscillations as a function of 
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FIG. 3: (Color online) Comparison between the energy split- 
ting of the Majorana mode at constant chemical potential 
(orange lines) and constant density (blue lines). The quasi- 
periodicity of the oscillations at constant \i is absent in a wire 
with fixed number of particles. The difference is particularly 
striking in the single-band case (top panels). Under the con- 
stant density condition, small variations of the initial occu- 
pancy result in qualitative changes of the dependence of AE 
on Vz- In contrast, the constant chemical potential condition 
is characterized by a generic oscillatory splitting. 

the magnetic field. 

The basis of our mechanism is that the Coulomb re- 
pulsion among the carriers in the nanowire, even after 
screening by the superconductor, renormalizes the elec- 
trostatic potential profile in the nanowire as a result of 
a self-consistent change in the density of electrons. Be- 
cause of the repulsive sign of the Coulomb interaction, 
this extra electrostatic potential will oppose any change 
of the electron density that might result from variations 
of the Zeeman potential Vz- Considering, for simplic- 
ity, the extreme limit of strong Coulomb interaction, we 
obtain the condition that the total number of electrons 
in the wire mjust be fixed, i.e. the nanowire behaves 
as a 'constant density' rather than a 'constant chemi- 
cal potential' system. In the numerical calculations we 
impose this condition by self-consistently adjusting the 
chemical potential so that the total number of electrons 
be the same as the particle number at the TQPT for a 
nominal non-interacting chemical potential A/iq. The 
results are shown in Fig. [3j The striking difference be- 
tween the dependence of the splitting on Vz at constant 
chemical potential versus constant density is particularly 
evident in the single-band case (top panels). To under- 
stand the suppression of the oscillations in the constant 
density case, we note that a single-band system in the 
topological phase has one occupied spin-channel. Con- 
sequently, the electron-density in this regime is tied to 
the Fermi wave — vector kp, e ff ~ n - Since kp, e ff dom- 
inates the rapidly oscillating part of AE in Eq. [3j if \i 
self-consistently adjusts to Vz in a way to keep kp, e ff 
constant, then AE becomes slowly varying in Vz as well. 

In the more realistic (and necessarily more complex) 
multi-band case, the particle density is not tied to the 
effective Fermi wave — vector of the top band, so one 
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FIG. 4: (Color online) Left panels: MF splitting as a function 
of chemical potential \i and Zeeman field Vz for single-band 
system (top) and for a wire with n top = 3 (bottom). The 
dark blue regions correspond to the minima of the energy 
splitting. The parallel bands represent the source of the oscil- 
latory splitting and a hallmark for Majorana physics. Right 
panel: Density as a function of chemical potential \i and Zee- 
man splitting Vz- Note that, in the single-band case, the con- 
stant density contours are nearly parallel with the structures 
in the left panel, which may results in a strongly suppressed 
splitting or in a monotonic dependence of the splitting on Vz- 

expects the oscillatory behavior to be partly restored. 
Nonetheless, the quasi-periodicity of the oscillatory split- 
ting is not a generic feature in a wire with constant par- 
ticle density. Specific examples are shown in the lower 
panels of Fig. [3] corresponding to a system with n top = 3 
(i.e., five occupied spin sub-bands). In systems with 
more occupied bands, the dependence of the splitting 
on V z approaches the quasi-periodic oscillatory behav- 
ior characteristic of the constant chemical potential case. 
To gain a deeper understanding of the manifestations of 
the Majorana oscillatory splitting in different conditions, 
we show in Fig. [4] (left panels) the dependence of AE 
on both the Zeeman field Vz and the chemical potential 
A//, together with contour plots of the electron number as 
function of the same variables (right panels) . Examining 
Fig. [4] we note that both the constant density contours 
and the characteristic parallel structures characterizing 
the MF splitting slope in a similar direction. Therefore, 
the MF splitting as a function of the Zeeman field at con- 
stant density will show far slower oscillations, than the 
constant chemical potential situation. Moreover, a con- 
stant density path may lie within either a low-splitting 
or a high-splitting band for a significant range of Zeeman 
fields, thus explaining the types of behavior illustrated in 
Fig. || 
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Our work clearly establishes that the Majorana- 
induced zero bias conductance peak should manifest 
some oscillatory splitting at large applied magnetic fields 
(or as a function of the gate-tuned chemical potential at 
constant magnetic field), and a direct observation of this 
oscillatory splitting behavior could be the smoking gun 
evidence for the Majorana existence. 

This work is supported by Microsoft Q, DARPA- 
QuEST and the Harvard Quantum Optics Center. 
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Numerical methods 

We consider a semiconductor nanowire with rectangu- 
lar cross section and dimensions L x ^> L y ~ L z . In the 
numerical calculations we have L y « 90nm, L z « 60nm, 
and two wire lengths, L x = l//m and L x = 2fim. For 
an InSb nanowire with this cross section, the typical en- 
ergy gap between the subbands is of the order 5mev. We 
assume that only a few subband are occupied and we 
project the single particle quantum mechanics problem 
for the nanowire electrons onto a low-energy subspace 
spanned by the states with energies lower than a certain 
threshold (typically ~ lOOmeV). As a result of the prox- 
imity to an s-wave superconductor (SC), a pair potential 
A is induced in the nanowire and the energy scale for the 
quantum states in the semiconductor (SM) are renormal- 
ized. To account for this proximity effect, we integrate 
out the SC degrees of freedom and incorporate them as 
a surface self-energy term of the form [29] 



E( W ) 



+ (T Z 



■ W 



(4) 



where 7 = 0.3meV is the effective SM-SC coupling, r x 
and r z are Pauli matrices in the Nambu space, Ao = 
1.5meV is the pair potential of the bulk SC, and C is a 
proximity-induced shift of the chemical potential. In the 
present calculations we take £ = 0. Within the static ap- 
proximation \/Aq — uj 2 — >• Aq, the self-energy becomes 
E(cj) ~ — 7CJ/A0 — J&yTy and the low-energy physics 
of the SM nanowire with proximity- induced SC can be 
described by an effective Bogoliubov-de Gennes Hamil- 
tonian. This approximation is valid, strictly speaking, 
at energies much lower than Aq, but represents a very 
good approximation even for E ~ Ao/2. Explicitly, the 
matrix elements of the effective BdG Hamiltonian can be 
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written as 

H BdG {n,n) = Z[e 

nn ' z 

+ Z{Hl ol ) nn t + A<7 y T y , (5) 

where n = (n xi n yi n z ) are quantum numbers for the 
nanowire states in the absence of spin-orbit coupling, e n 
are the corresponding energies, Vz is the Zeeman split- 
ting, and (H§ OI ) nn f and (Hg 01 ) nn > are matrix elements 
for the intra-band and inter-band Rashba spin-orbit 
coupling, respectively. Note that energy scale for the SM 
nanowire is renormalized by a factor Z=(l + 7/Ao) _1 
due to the SC proximity effect. This renormalization is 
determined by the term in the self-energy Q that is pro- 
portional to uo (in the static approximation). The pairing 
term in Eq. ([5| is derived from the corresponding con- 
tribution to the self-energy Q and is proportional to the 
induced pair potential A = 7Ao/(7 + Ao) = 250/ieV. 
The effective Hamiltonian described by Eq. ([5| is diago- 
nalized numerically. 

Analytic solution for Majorana wave-functions in 
nanowires 

The one-dimensional BdG Hamiltonian for a single 
band semiconductor nanowire with spin-orbit coupling 
(assumed to be linear in the momentum k y ), which in 
general can be written as 

H B dG = (~dl - fi(x))r z + V z a z + iad x a y r z + Ar x (6) 

where Vz is the strength of the effective Zeeman field and 
the unit vector a characterizes the spin-orbit coupling. 

Non-degenerate Majorana spinor solutions are of the 
form \I/ = (u, idyll*) and are completely determined by 
the 2-spinor u. This fact was used to obtain the Majo- 
rana solutions for vortices to reduce the BdG equation 
from a 4 x 4 system of equations to a 2 x 2 system of 
equation. However, this reduction procedure required the 
BdG Hamiltonian to be real which is not the case for gen- 
eral forms of spin-orbit coupling and Zeeman splitting. 
The BdG equation for the zero energy mode HsdG^ — 
may be reduced to an equation for u as 

[(-dl - n(x)) + V Z (J Z + a(ia x )d x ] u + A(ia y )u* = 0. 

(7) 

This equation is not real but may be reduced to a system 
of real equations by writing u = ur + iuj and taking the 
real and imaginary parts of the resulting equation giving 
a pair of equations of the form 

[{-dl - n(x)) + V Z (Jz + a(ia y )d x + A(ia y )} u R = 

(8) 

[(-dl - fjL(y)) + V z <Tz + ot(ia y )d x - A(ia y )] m = 0. 

(9) 



This results in a reduced BdG equation for the E = 
reduced spinor ^(x) 

\ -AA - ad x -dl -V z - fi(x) J v 7 v 7 

where A = ±1. 

An end of the wire of the type considered above is 
defined by requiring some parameter of the Hamiltonian 
to vary across the end situated at x = 0. We take this 
parameter to be constant for x < and x > 0. In this 
case, our previous approach can be applied in a way even 
simpler than the application to the vortex problem, since 
the solutions on both sides of the interface at x = can 
be approximated as a sum 

U ( X ) = ^2 a ne~ ZnX pn (11) 

n 

where 

( ~zl + V z - n XA-z n a \ . 

V -\A + z n a -zl-V z -n) Pn -°- (12) 

The solutions z n must satisfy the quartic equation 

Det(- Z \ + A V : -» X 2 A ~ Za ) 
\ -AA + za -z A —V z —fi J 

= (z 2 + - V z 2 + (za T A) 2 = 0. (13) 

For 

C = (A 2 + M 2 -y|)<0, (14) 

there are 3 values of z n such that Re(z n ) < in a given 
A channel. The coefficients a n in the solution are deter- 
mined by matching the boundary conditions on ^(x) at 
x = 0. 

The equation Eq. [l3j while in principle has an exact 
solution, the solution is extremely complicated to write 
out explicitly. Therefore, we will discuss an approximate 
solution, which will be valid in the limit where either a 
or A are small. In this case Eq. [13] is written as 

(z 2 + fi) 2 - V 2 + (za =F A) 2 
= (z 2 + /i) 2 - V 2 + (z 2 a 2 =F 2aAz + A 2 ) = 
= z A + (a 2 + 2fi)z 2 - (V 2 - A 2 - ii 2 ± 2aAz) = 

z = ± 

\J-(H + a 2 / 2) ± y/(fjL + a 2 /2) 2 + (V z 2 - A 2 - /i 2 ± 2a Az). 

(15) 

In the limit where spin-orbit coupling and a or A is 
small, the gap in the bulk BdG spectrum is small and 
one can have nearly propagating solutions, which are ap- 
proximated as below. Since the values of z determine the 
x-dependence of the solutions of Eq. [TT] z can be written 
as 

z = ±ik Fie ff - £~\ (16) 
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where £ is the effective coherence length and kp, e ff is the 
effective fermi- wave- vector associated with the zero-mode 
solution *ff(x) oc e lkF > e ff x ~ x /€ for x > £. In the limit 
where either of a or A are small, £ can be approximated 

by 



r 1 



2y/{ti + a 2 /2) 2 + Wz " A2 " M 2 ) 
and kF, e ff can be approximated by 



(17) 



k F,eff 



'z-A 2 -M 2 ) + (M + « 2 /2). 

(18) 

To determine the full Major ana wave- function in 
Eq. [TT] in addition to z n one needs to determine p n = 
(i£ n ,t> The spinors p n are calculated by solving the 

equation 



-z 2 + V z -p 
-AA + za 



-z 



AA — za 

2 -v z - 



The solutions are now given by 

Uf = — (AA — za) 



V z 



/i, 



= 0. (19) 



(20) 
(21) 



where z is an acceptable solution for the quartic equation 
Eq.pl (Eq.pl, for V 2 > (A 2 +p 2 ) there are 3 solutions 



on the right half of the complex z plane and 1 solution 
on the left half for A = —1. The situation is opposite for 
A = 1. Since the equations are real, the 3 solutions are 
z, z* and w where w is real. The wave-function ip(x) is 
then written as 



u{x) = e 



(AA - pa) 

2 ^V 2 



a e 



-(AA - z*a) 
z* 2 + V, 



p 



ae 



p 



-(AA - za) 
-z 2 + V z - p 

(22) 



The boundary condition ^(0) = implies that 

_ i(z* - w)[V z a - z*(\A - aw) - (an + AXw)] 



2Im(z)(V z a + \z\ 2 a - 2Re[z]A\ - an) 



(23) 



Since the four components of the spinor are related by 
particle- hole symmetry, ip = (u, ia y u*), we can write i(j 
as 



The left Majorana fermions ^i(x) satisfies 



[{-d 2 x -p- id x cr y )r z + V z a z + Ar x ]^i(x) = 0. (25) 



An inverse solution is obtained by noting that ip r ( x ) = 
ia z ipi(—x) is a solution as well, which decays in the op- 
posite direction. The factor of i ensures that ijj r (x) is 
invariant under a y r y K. 



Splitting of MFs in spin-orbit coupled nanowires 



To compute Majorana fermion splitting we proceed us- 
ing the variational principle by writing 



V(x) = ciijji(x) + c r ijj r (x), 



(26) 



where ip^ r (x) are the Majorana fermion wave-functions, 
which are zero-modes of the Hamiltonian for x > and 
x < L respectively that have been appropriately trun- 
cated at the other boundary. Here Sk(x) is expected to 
be a solution for a wire with chemical poetntial domain 
wall at both ends x = 0, L. 



E 



j dx^(x)^(x) 
c*c r J dxtl)l(x)Htp r (x) +c*ci J dxtl)*(x)Htpi(x) 
|q| 2 + M 2 + c z *c r J dxil)l(x)t/j r (x) + c*q J dxip*(x)il>i(x) 

c\c r I dxip* (x)Hip r (x) + c*q / dxip*(x)Hipi(x) 



nL/2 pi 

\c r \ dxipl(x)Hi/j r (x) + c*ci / dxip*(x)Hipi(x) 
Jo Jo 

L 
L/2 



r L/2 



p Lj n Lt 

+ c*c r / dxipi (x)Hifj r (x) + c*ci / dxip*{x)Hipi{x)^ 

J L/2 ' ' J L/2 

(27) 



ip(x) — e 



( -(AA - wa) \ 
-w 2 + V z - p 
-w 2 + V Z — p 
\ (AA - wa) J 

( -(\A-z*a) \ 
z* 2 + V z - p 
z* 2 + V z -p 
\ (\A-z*a) J 



ae 



( -(AA - za) \ 
-z 2 + V z - p 
-z 2 + V z - p 

\ (AA - za) J 

(24) 



where we have assumed |q| 2 + \c r \ 2 — 1 and also used 
the fact that ipi jr (x) are particle- hole symmetric (so 
their expectation value for H vanishes). Observing that 
Hi/j r (x) = for x > L/2 and H^ t (x) = for x < L/2, 



pL/2 r>iu 

iC r / dxijji(x)Hip r (x) + c*ci / dxip*(x)Hipi(x). 

Jo Jl/2 

(28) 



E w c] 
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Using Green theorem, where K(x — x') is the metallic screened Coulomb kernel 

b b and p(x f ) is the density of electrons on the nanowire. 

/ dx^{x)d x (j){x) = / dxd x {^{x)(j){x)) - d x $* {x)<j>{x) Within BCS theory, the total charge density p(x') is given 

J a J a \yy 



[(fTO)la- / dxd x r(x)^2 
J a 



(29) 



p(x) 



^2 \u n (x,a)\ 2 , 



(35) 



dxip* (x)d 2 (j)(x) = 
W{x)d^(x))] b a - 

b 

dxd 2 x ^(x)(j){x) 



n:E n <0,a 



dxd x {^*{x)d x (j)(x))-d x ^*(x)d x (i){x\ . 

where n is summed over states with h n < 0. As a consis- 



dxd x {d x ijj* (x)(j)(x)) 



= [{r{x)d x ^x))-d x r{x)4>{x)t+ f dxdirwm 

(30) 

L/2 

dxi/j*(x)[(-d% - p + iaa y d x )r z + V z a z + Arjf (j) 

r L/2 



/ dx[{(-c^ - /i + iaa y d x )r z + + Ar,}^(x)]^(x] 
Jo where ft 



tency check note that in the absence of superconductivity 
J dx^2 a \u n (x,a)\ 2 = 1 for states with (e n — /i) < and 
vanishes for other states. Therefore J dxp(x) is the total 
number of electrons for A = 0. 

For simplicity, we ignore the spatial variation of p(x) 
and approximate 

p(x) « p (36) 
and also ignore the variation in V ex t(x) so that 

Vext{x) « Kxt = « _1 P, (37) 



is the mean compressibility. 



+ [(^(x)d x r z i/; r (x) + iai/ji (x)a y T z ilj r (x))}Q /2 

= (^(L/2)d x r z ML/2)^ia^(L/2)a y r z ML/2)). 

(31) 

Using the relation tp r (x) = a z tpi(L — x), we note that the 
above matrix element becomes 

(tf(L/2)d x T z a z ML/2) + iari{L/2)a v T z cj z iPi{L/2)). 

(32) 

In the case of a long wire in the case where one of the 
modes has a significantly longer decay length 



The mean-external potential V ex t renormalizes p in the 
BdG equation according to 



p -> p = p - V ec 



(38) 



Assuming that the compressibility is small ft — » 0, it is 
clear that a small change of p changes V ex t and therefore 
mu by a large amount. This leads to the constraint that 
when parameters other than V ga te are changed 



» J dxp(x) 



0. 



(39) 



E^h 2 k F - 



■ cos (kpL), 



(33) 



where kp^ are given by Eq. 17 and Eq. 18 This equa- 
tion appears as Eq. [3] in the main manuscript and is the 
central result of this section of the Supplementary mate- 
rial. Based on Eq. [T8j kp depends on Vz, which enters 
into a rapidly varying (in the limit of large L) cosine term 
in Eq. [33] Therefore, one can expect the splitting E to 
oscillate as a function of Vz as discussed in the main text. 



It is possible to change p(x) by changing the gate voltage 
- but as in mesoscopic wires, one needs to change V ga te by 
many volts before a change in density corresponding to 
a few meV change of chemical potential is accomplished. 

In the simple single-band case and ignoring supercon- 
ductivity A, the linear density of electrons is given by 
n = kp/2 and the chemical potential p is given by 



p = 72,74- JV 2 -n 2 a 2 /4. 



(40) 



Role of Coulomb interactions 
calculations 



self-consistent 



Substituting this expression for p into Eq. 18 we find that 



Let us now suppose that the electrons in the wire are 
subject to a Coulomb interaction that is screened by the 
superconducting substrate by perfect metallic screening. 
Within the Hartree (the Fock term is ignored here) ap- 
proximation, one expects a renormalization of the exter- 
nal potential 



V ex t(x) 



J d x>K ia 



x')p(x'), 



(34) 



in the limit A — > 0, kp y eff ~ n/2, which is independent 
of Vz- Since based on Eq. [3J the oscillations as a func- 
tion of Vz are a result of the Vz dependence of kp, e ff, 
the oscillations in the Vz dependence can be expected 
to be suppressed in the constant density limit. This is 
confirmed in the single-band case by results in Fig. [I] 

The multi-band results are more complicated and in 
principle can depend on the inter and intra band capac- 
itances of various channels. However, assuming that all 
the capacitances are the same, the multi-band effect is 
expected to screen the role of the coulomb interaction, 
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so that changing V z leads to some change of hp. The The Vz dependence of the splitting in the multi-band 
result (shown in Fig. [T]) leads to more ^-dependence of constant density case can be understood by examining 
the splitting than in the single-band case. Fig. [2] 



